**FULL TITLE** 

ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** 
**NAMES OF EDITORS** 



The small-scale solar surface dynamo 



Jonathan Pietarila Graham, Sanja Danilovic, and Manfred Schiissler 

M ax- Planck- Institut fiir Sonnensystemforschung, 37191 
Katlenhurg-Lindau, Germany 

Abstract. The existence of a turbulent small-scale solar surface dynamo 
is likely, considering existing numerical and laboratory experiments, as well as 
comparisons of a small-scale dynamo in MURaM simulations with Hinode obser- 
vations. We find the observed peaked probability distribution function (PDF) 
from Stokes-y magnctograms is consistent with a monotonic PDF of the ac- 
tual vertical field strength. The cancellation function of the vertical flux density 
from a Hinode SP observation is found to follow a self-similar power law over two 
decades in length scales down to the w 200 km resolution limit. This provides 
observational evidence that the scales of magnetic structuring in the photosphere 
extend at least down to 20 km. From the power law, we determine a lower bound 
for the true quiet-Sun mean vertical unsigned flux density of w 43 G, consistent 
with our numerically-based estimates that 80% or more of the vertical unsigned 
flux should be invisible to Stokes— y observations at a resolution of 200 km ow- 
ing to cancellation. Our estimates significantly reduce the order-of-magnitude 
discrepancy between Zeeman- and Hanle-based estimates. 



1. Introduction 

The existence of a global solar dynamo is generally accepted, though it is not 
well understood. There is evidence for a second, loc al dynamo operating near t he 



solar surface - the so-called solar-surface dynamo (jPetrovav Sz Szakalvlll993l )r 
This dynamo may be driven by turbulence: magnetic energ y is amplified b y 
random stretching of magnetic field lines by turbulent motions ( Batchelor 19501 ) 



Magnetic energy is lost to Ohmic dissipation. For dynamo action to succeed, 
stretching must dominate over dissipation. The relative strengths of these two 
effects is quantified by the magnetic Reynolds number, Rcm = vqIq/t] [vq and Iq 
are typical velocity and length scales and r] the magnetic diffusivity) . Similarly, 
the kinetic Reynolds number is Re = vqIq/v where v is the kinematic viscosity. 
When ReM exceeds a critical threshold, ReM > Re-%[, dynamo action results. 

Two different types of dynamo can be found, depending on the presence or 
absence of net flow helicity (jMeneguzzi et al.lll98ll ). With net helicity, magnetic 



energy grows at scales larger than the energy-containing scale of the fluid mo- 
tions: large-scale dynamo (LSD) or mean-field dynamo. LSDs are often studied 
with mean-field theory; the production of large-scale magnetic energy is approx- 



^ Actually, in simulations, turbulent dynamo action occurs in the bulk of the convection zone, 
though greater t han solar rotation rate may be required to generate large-scale field in the 
absence of shear (|Brown et alj[2007h . 
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imately the alpha-effect (see, e.g., iBrandenburd [20031 ) . Without net hehcity, 
dynamo action is harder to achieve and magnetic every grows at scales smaller 
than the forcing scale. This latter defines small-scale dynamo (SSD) or fluctu- 
ation dynamo action. Near the solar surface, the convective time scale is much 
shorter than the rotation period, the effects of rotation can be neglected, and a 
flow with no net helicity results. Any surface dynamo will thus be a SSD. 

This is also suggested by observation of the small-scale magnetic field in 
the quiet Sun. In high resolution magnetograms we see mixed polarity fields on 
small sca .les, which is variously called the magnetic car pet or the salt-and-pepper 
pattern ( Title &: Schriivei]|l998l : [Hagenaar et al.|[2003 ). consistent with the idea 
of SSD action. As turbulent convection can drive small-scale dy namo action (in 
numerical simulation s of Boussinesq convection without rotation ICattaneolll999l : 
Cattaneo et al.|[2003l ) , observations and simulations together provide evidence of 



a SSD driven by turbulence at the solar surface (likely deeper as well). 

Several arguments can be put forward against a small-scale solar dynamo. 
Firstly, there exists the possibility that the small-scale field is solely produced by 
the shredding up of large-scale field by turbulence. However, ob servationally the 
amount of small-scale fiux is no t dependent on the solar cycle ( Hagenaar et al.l 
20031 : iTruiillo Bueno et"aDl2004l ^. This might not be the case if the small-scale 



flux is the result of the shredding up of the field from the global dynamo. It might 
also then show some latitudinal dependence (among low latitudes). Assuming 
the existence of both SSD and shredding, small-scale dynamo is predicted to 

1 /2 

create sn iall-scale magnetic field at the turbulent rate of stretching (a Rej(,j , 
see, e.g., Ilskakov et al.l 20071 ) which is much faster either than large-scale field 
can be produced (at time scales associated with t he kinetic-energy-containi ng 
length scales) or be shredded up by the turbulence ( Schekochihin et al1l2005l ). 

Secondly, small-scale dynamo action may not be possible at the magnetic 
Prandtl number, Pm, of the solar plasma, Pm = RcM/Re ~ 10~^: -Re^j 
sharply increases with decreasing Pm (increasing Re) since eddies smaller than 
the characteristic scale of the magnetic field diffuse the field and inhibit dy- 
namo action. T he two asymptotic possibilities are Re^^ — >■ const as Re — > ■ oo 
(SSD at low PA^: lRogachevskii k Kleeorinlll997l : iBoldvrev k Cattanedl2004ll or 



Schekochihin et al 



Re^/Re = P^ ^ const as i?e — > oo (no SSD at low Pm 
I2OO5I ). Numerical simulations of low Pm SSDs have focused on the existence 
(or not) of a time-averaged mean flow. Su ch a flow exists f or the Sun and 
many other astrophysical cases . Both with ( Pontv et al. 20051 ) and without a 
mean flow ( Iskakov et al.l [20071 ) . a plateau in Re'^ was found. This suggests 
such dynamo action should be possible on the Sun. However, the Pm of the 
sun is 3 orders of magnitude smaller that that accessible to present numerical 
computations; the observed plateaus may not represent the asymptotic behav- 
ior. Fortunately, a laboratory dynamo resulting from u nconstrained turbulence 
in liquid sodium {Pm ~ 10~^) has been demonstrated ( Monchaux et al.ll2007l ) 
establishing that a turbulent dynamo is possible at values of Pm corresponding 
to the solar plasma. 
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2. The MURaM code and the solar surface dynamo 



The r emaining; o bjection against a surface dynamo was raised by IStein et"al 



( 20031 ) : unlike the lCattaneo ( 19991 ) simulation, the Sun is strongly stratified and 
magnetic flux is swept into the d own-flow lanes an d subject to long recircula- 
tion times. In the simulations of lStein et al. ( 20031') with open b oundaries, no 
dynamo action was found. However, Vogler &: Schiissleil ( 20071 ) have demon- 
strat e d SSD action for s ufficiently high Rem using the MURaM code (jVogler 
l2003l : IVoder et al.ll2005l l. They simulated a box with strong stratification and 
open boundaries, but prevented advection of magnetic fiux from outside the 
box: providing an artificially isolated surface layer ignoring any SSD action 
that may occur in deeper layers. They conducted 3 simulations with increasing 
Rcm (see Tabled]). Above the critical threshold, Re'^j ~ 2000, dynamo action 
occurs. The magnetic energy spectrum peaks at scales smaller than the energy- 
containing scale of the fluid motions, demonstrating that SSD is a possibility 
for a solar surface dynamo despite strong stratification and little recirculation. 



In Ta ble [H dynamo action in the 3 MURaM si i nulati ons of lVogler &: S chiissleil 
(120071 1 as well as in the simulations of lCattaneol (j 19991 ) and of I Stein et al.i ((200j) 
is indicated solely by their Rcm and Pm- For the simulations below the crit- 
ical threshold for dynamo act ion, {Re^ ^ 1000), there is no dynamo action. 
The result by ICattaned (|l999l ') is an exception as Pm > 1 and this lowers the 
threshold. 



Table 1. Simula t ions by ICattaneol (|l999[) . IStein et al.l (|2003l ). and 



Vogler fc Schiissleil (|2007l) : computational grid size, magnetic Reynolds and 



Prandtl nmnbers, boundary conditions (BC; open or closed), and presence of 
SSD. 



Simulation 


Grid pts. 


Reu 


Pm 


BC 


SSD 


Run A 


288^ X 100 


« 300 


< 1 


open 


N 


Stein et al. 


253^ X 163 


« 600 


w 1 


open 


N 


Run B 


576^ X 100 


« 1300 


^1 


open 


N 


Cattanco 


5122 X 97 


« 1000 


w 5 


closed 


Y 


Run C 


648^ X 140 


« 2600 




open 


Y 



The MURaM Run C dyn amo also reproduc es, quantitatively, some aspects 
of observations. For instance, iLites et al.l ( 20081 ) find from Hinode SP observa- 
tions that the mean horizontal magnetic field is stronger than the mean unsigned 
vertical magnetic field. Using a spatial map, in the "normal mode" (324 xl64 ), 
they make magnetograms of the apparent longitudinal flux density, i^app' ^^'^ 
apparent transverse flux density, i?app- They flnd a factor of 5 stronger hori- 
zontal fields: the means are (\B^:r^T^\) « IIG and ( 55G. A similar 



ratio w as found in the MUR aM Run C dynamo bv ISchiissler Sz Vogler ( 20081 ) : 
see also Steiner et al.l ( 20081 ) . In the MURaM simulation, both small-scale, low- 
lying vertical loops and extended canopy-like structures are found in the line 
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formation region. These structures (related to flux expulsion of horizontal field) 
contribute to the stronger mean horizontal fields: horizontal field occupies a 
larger area than the narrow vertical foot points. The simulation also makes a 
prediction: lines with different formation heights will yield different ratios. 



3. Measurement of the turbulent magnetic field 



Whatever its source, the small-scale quiet-Sun magnetic field is t urbulent. This 

shou l d be taken into consideration whe n interpreting observations ([Sanchez Almeida et al, 
Il996l : ISanchez Almeida fc Lited I2OO0I ) . Consider the discrepancy between the 
probability distribution functions (PDFs) of the apparent vertical flux densities 
from the Hinode SP observation (discussed in the previous section) and that of 
B;, in MURaM dynamo Run C (see Fig. [I]). For the SP observation, the PDF 
is peaked at ~ 3 G. The simulation PDF is monotonic; there is a much greater 
amount of weak field than in observations. Here -Bave is defined as averaged 
over the height range corresponding to logr € [—3.5, .1]: other samplings show 
similar PDFs. Synthetic Fel 630 nm doublet spectra are generated from the 
MURaM cube and used to create magnetogram signals. The PDF of the result 
is shown as a dashed line in Fig. [T]and begins to resemble the peaked PDF of the 
SP observation. With the addition of instrumental noise the PDF (dash-dotted 
line) becomes quite similar to that of the Hinode observation. The observa- 
tion has a much greater amount of strong field due to the much higher Rcm of 
the Sun (more efficient SSD) and the presence of network elements which are 
not in the simulation. In any case, these considerations show that the observed 
peaked PDF is in fact consistent with a much greater prevalence of weak vertical 
magnetic field. 




(G) 



Figure 1. PDFs for magnetic field strengths and derived proxies: MURaM 
simulation Save (solid line), Hinode SP observation i?app (dotted line), and 
synthetic i?app from MURaM (dashed) - including noise (dot-dashed). 
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Vertical radiative t ransfer through a turbulent fl uid contributes to this de- 
formation of the PDF ( Sanchez Almeida et al.ll2003l ). In Fig. ^a), we plot the 
z— components of the magnetic field and velocity for a simulation pixel with 
Bave ^ -^app- The magnetic field undergoes a reversal above r = 1. The 
average of over the line formation region is very smah (~ 10"^ G), but the 
reversal of Bz is coincident with a strong gradient in v^- Thus, the absorption 
profiles from the two opposite-polarity fields are Doppler shifted with respect 
to each other and do not cancel. Such cancellation leads to a stronger Stokes 
V signal than would otherwise arise (also an asymmetric one, see Fig. E^b)). 
If this statistically accounts for the differences in the PDFs, a correspondence 
between the strength of the velocity fluctuations along the line of sight and the 
synthetic magnetogram would be expected (for pixels with weak i^ave)- This 
is indeed seen in Fig. [3l which includes only pixels with |-Bave| < 0.1 G. There 
exists a definite trend of stronger magnetogram signal when there are stronger 
Doppler shifts (velocity gradients) between the various layers of the atmosphere. 




Figure 2. (a) Bz (solid line) and Vz (dashed line) versus rsoonm, and (b) 
Stokes V profile for a pixel with -B^pp = — 9G and Save = ^5 • 
Strong gradients can lead to |-Bapp| ^ | -Save I- 

How do gradients in the horizontal direction affect observations? Or, given 
that opposite-polarity vertical field inside a resolution element leads to cancella- 
tion of Stokes V signal, how can the turbulent (and hence fractal) nature of the 
magnetic field be used to estimate the cancellation (i.e., to get an estimate for 
the true unsigned vertical fiux density)? Fractals are self-similar; they display 
power-law scaling such as in the box-counting fractal dimension: the domain is 
partitioned into boxes of edge length I and the number of boxes containing the 
fractal, N{1), has a power-law relation, N{1) oc , where Df is the fractal 

dimension. The magnetic field, however, is more complex; for each box in a 
partition, there is a net flux (a magnitude) and a direction (sign), either up or 
down. In this case, we need to use the cancellation function, 



x{i) ^ E 



/ Bzda 


/L 







\Bz\da (1) 
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Figure 3. Average -Bapp versus the fluctuations of tire vertical velocity 
along the line-of-sight, cry (|i?ave| < 0.1 G for all pixels). 



introduced by lOtt et al.l ( 19921 ) . Simply put, measures the portion of net 
flux remaining after averaging over boxes of length /. A fractal (self-similar) 
magnetic field has a power-law scaling, x(0 ~ -^or the magnetogram from 
the Hinode SP observation (see Fig. H] (a)), a clear power-law scaling of the 
cancellation function exists over two decades of length scales down to the reso- 
lution limit at ~ 200 km. Due to experience with turbulent power-laws, we can 
be confident that more cancellation occurs at scales do wn to ~ 20 km or less as 
the power-law continues down to the resolution limit ( Carbone &: Brunol 1997 



Sorriso-Valvo et al.l[2003 ). It is unknown how far this power law extends, but it 



must stop before the magnetic dissipation scale, Assuming that the power 
law continues until dissipation sets in, we derive an estimate for the true, mean 
unsigned vertical magnetic field for the Sun, {\Bz\), defined as 



(l^.l) 



da . 



The mean absolute value of measured at resolution / is 



/ Bzda 




JA^{l) 





(2) 



(3) 



Note that = x(0 ' (l^zD- Using that there is no cancellation below the 

magnetic dissipation scale, {\Bz\)i^ = (j-B^I), 



{m) = {\Bz\)i 



x{i) 



12G- 



100 km \ 



0.26 



(4) 



where 12 G is the 100 km result from lLites et al.l (j2008l ). We a pply Kolmogorov 



theor y to estimate for the Sun. Assuming rj 108 cm s i\ 



Kovitva &: Cram 



19831 ). we find Irj to be ~ 80 m. At this scale, x(/) must have a slope of zero and 
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the power-law is affected for about a decade of scales larger. There will still be 
some cancellation but the slope will be decreasing. Therefore, our calculation is 
based on a scale of 800m in Eq. giving a lower bound of ~ 43 G. 

An estimation for (l^^l) can also be obtained using the MURaM dynamo 
simulation. In Fig. Il](b), we plot the cancellation at a resolution of 200 km ver- 
sus Rsm. The power-law scaling can be extrapolated to the expected Rcm of the 
Sun ~ 3 ■ 10^: we estimate x(200A;m) ^1/5, which means (ji^zl) ^ 11-5 ~ 55 G. 
This provides qualitative agreement with Hanle estimations. Based on our esti- 
m a.tes, the mean uns igned vertical field strength is about 40 or 50 G. As shown 
bv lLites et all (|2008l 1. the mean horizontal field strength is stronger. Therefore, 
the mean vector magnitud e of the magnetic field of ~ 130 G as reported by 
iTruiillo Bueno et al.1 ((2003 ) is not in contradiction with these Zeeman results. 




Figure 4. (a) Normalized cancellation function, x{l) /x(l Mm) , versus scale, 
I, from Hinodc i?app- The exponent of the fitted power law is k = 0.26±0.01. 
(b) Portion of flux remaining at / = 200km, x(200km), versus Rcm for 
MURaM dynamo. A power-law is indicated by the dashed line. 



4. Summary 

A small-scale solar dynamo near the surface is likely. MURaM simulations show 
dynamo action and its properties are found to be in agreement with Hinode 
observations. The arguments against a solar small-scale dynamo, thus far, have 
failed. Whatever its source, the small-scale magnetic field is turbulent and fractal 
and this should be taken into consideration when interpreting observations: 1) 
the PDF of Stokes V field estimates do not accurately represent the PDF of the 
Bz and are consistent with a much greater prevalence of weak B^; 2) the multi- 
fractal self-similar pattern of the quiet-Sun photospheric magnetic field extends 
down to the resolution limit, 200 km. This constitutes observational evidence 
that the smallest scale of magnetic structuring in the photosphere is 20 km or 
smaller. The power law also constrains the quiet-Sun true mean unsigned vertical 
flux density: the lower bound, ~ 43 G, is consistent with estimates based solely 
on numerical simulations (~ 55 G). The order of magnitude disparity between 
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Hanle and Zeeman-based estimates may be resolved by a proper consideration 
of the cancellation properties of the full vector field. 
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